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Abstract 

We discuss stochastic representations of advection diffusion equa- 
tions with variable diffusivity, stochastic integrals of motion and gen- 
eralized relative entropies. 
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1 Introduction 

Recently, Michel, Mischler and Perthame jH] discovered a remarkable prop- 
erty of certain unstable linear equations, in which decay of relative entropies 
takes place. Their observation was applied to population dynamics models, 
but the list of applications is growing. Of course, relative entropies have 
been used for a long time in kinetic theory and conservation laws. However, 
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the decay of relative entropies, was known before only in stable, self-adjoint 
situations in which a global attracting steady solution exists and no flow 
advection is present jU]. The property of decay of relative entropies was 
slightly generalized to variable diffusion coefficients and applied to Smolu- 
chowski systems in [1.. A stochastic interpretation and proof in the case of 
constant diffusion coefficients was given in Here we provide a stochastic 
interpretation and proof in the case of variable diffusion coefficients. The 
method of proof and concepts are of more general interest [31 Hj. 
We consider a linear operator 

Vp = vdi(aijdjp) - div x (Up) + Vp (1) 

in lR n , where 

U{x,t) = {U j {x,t)) j=1 ^ n (2) 

is a smooth (C 2 ) function, V = V(x,t) is a continuous and bounded scalar 
potential and 

Oij(x, t) = cr ip (x, t)<T jp (x, t) (3) 

with the matrix 

a(x,t) = {o-ij(x,t))ij (4) 

a given smooth (C 2 ) matrix. We assume that a is bounded and U and V x a 
decay at infinity. We use the shorthand notation A(D) for the operator 

A(D)p = aijdidjp (5) 

and use also the non- divergence form 

Vp = vA(D)p - u ■ V x p + Pp (6) 

where 

Uj(x, t) = Uj(x, t) — i>di(aij(x, t)) (7) 

and 

P = V-div x (U). (8) 
The formal adjoint of the operator T> in L 2 (M. n ) is 

V*(j> = vd i {a ij d j <j>) + U- V x <j) + V<j). (9) 

The following is the result of Michel, Mischler and Perthame: 
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f H (-)(f)pdx<0. (13) 



Theorem 1 Let f be a solution of 

dtf = Vf (10) 
and let p > be a positive solution of the same equation, 

d tP = Vp. (11) 

Let H be a smooth convex function of one variable and let cp be a non-negative 
function obeying pointwise 

d t (j) + V*cf) = 0. (12) 

Then 

dt J \p 

2 Stochastic Lagrangian Flow 

In order to represent solutions of equations like (jlUj) we consider the drift 

Vj(x, t) = Uj + 2u(d k a jp )a kp = Uj - v{d k a kp )a }p + v{d k a jp )a kp . (14) 
Let X(a, t) be the strong solution of the stochastic differential system 

dXj(t) = Vj(X, t)dt + V2^a jp (X, t)dW p (15) 

with initial data 

X(a,0) = a. (16) 

Here is a standard Brownian process in M n starting at time zero from the 
origin. This process will be fixed throughout the paper and all measurability 
issues will be with respect to the filtration associated to it and all almost sure 
statements will be with respect to the probability measure on the standard 
Wiener space. We will need the following result: 

Theorem 2 The inverse of the flow map a X(a,t), the stochastic map 

x^A(x,t) (17) 
exists almost surely and satisfies its defining relations 

X(A(x,t),t) = x, Vx G 1", A(X(a,t),t) = a, Va G 1", Vt, a.s. 
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The map X is smooth and the determinant 

D(a,t) =det(d a X(a,t)) (18) 

obeys the SDE 

d(det (d a X(a,t)) = [det(d a X(a,t))]x 
U(div x v)(x,t) + 2vE(x,t)]\ x=x{at) dt + y/2v(d k (a kp ))(x, t)\ x=X ( a ,t)dW p } 

(19) 

with 

E(x, t) = J2J2 det (^k- (2°) 

i<j p 

The map A(x,t) satisfies the stochastic partial differential system 

dAj + (u ■ V x Aj - uA(D)Aj) dt + \f2v(d k Aj)a kp dW p = (21) 
with initial data 

A(x,0) = 0. 

Remark. In the statement above, det(diO-j P )ij refers to the determinant 
of the two-by-two matrix (d r a kp ) with r,k e {i,j} for fixed i < j and p. 
Theorem El was originally proved in [3] for constant coefficients and in jl] for 
variable coefficients. For completeness, we reproduce the proof (with variable 
coefficients, as stated above) in Appendix 1X1 

3 Stochastically Passive Scalars and 
Feynman-Kac Formula 

We consider first deterministic smooth time-independent functions fo and 
note that the functions 9 = 9f (x,t) = fo(A(x,t)) are stochastically passive 
in the sense that they obey the equation 

d6+(u- V x 9 - vA(D)8) dt + V2^d k 9a kp dW p = (22) 

with initial data 

6(x,0) = f(x). (23) 

Solutions of the SPDE form an algebra; in particular, products of so- 
lutions are solutions, a nontrivial fact due to the presence of the stochastic 
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term. The expected values of these scalars obey advection-diffusion equa- 
tions and do not form an algebra in general, if v > 0. We consider now the 
function 

I(a, t) = exp j jf P(X(a, s), s)ds J (24) 
where P(x, t) is given in (JHJ) and consider the function 

i) = if> fo (x,t) = % (x,t)/(A(x,t),t) (25) 

We have 

Theorem 3 The process ip — ipf given by 

i/>(x,t) = f (A(x,t))exp{ I P{X{a,s),s)ds\ 

a=A(x,t) 

(26) 



solves 

dip - (Vip) dt + V2uV x ipadW = (27) 
with initial datum ip(x,0) = fo{x). 

The proof of this result follows using stochastic calculus [Sj, [Sj. Indeed, the 
function I(a, t) obeys 

d t I(a, t) = P{X{a, t),t)I(a, t) (28) 

pathwise (almost surely). Then, a calculation using (j2~I|) (see [I]) shows 
that the function 

J(x,t) = I(A(x,t),t) (29) 

solves 

dJ + (u ■ V X J -PJ - vA{D)J)dt + v^V, JadW = 0. (30) 
The function ipf is the product 

and therefore, from Ito's formula 

c% = Jde + 0dJ + d(J,9) 
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and the equations obeyed by J, 6, we have 

<%> = 

(-u ■ V x ip fo + Pip fo + uJA(D)9 + u9A(D)J + 2v(d k J)a kp (d j 6)a jp )dt 

-V^v x ip fo <Tdw. 

This means 

dip = (-u ■ V x i> + Pip + vA(D)ip)dt - V2vV ' x ip fo a dW. 
Because of (0) we have (J2*7)l . 

4 Stochastic Integrals of Motion. 

Proposition 1 Consider a deterministic function cf) that solves \lWt). Then 
the junction 

M{a,t) = <p{X(a,t),t)det(d a X(a,t))exp |jf P(X(a, s), s)ds| (31) 

zs a martingale. 

Proof. We start by writing 

M(a,t) = $(a,t)I(a,t)D(a,t) 

with 

$(a,t) = (j)(X(a,t),t), 

I given above in (|24)l and D given in p8jl . Next, we compute the equation 
obeyed by $/. In view of (J2*8j) and using Ito's formula we have 



d($I) = I {{d t <f>{X{a, t),t) + P(/)(X(a, t),t)) dt + 

+ V x <f)\ X (a,t) ■ dX + \didj(j)\x{a,t)d{Xi, Xj)} , 

which gives, in view of (fT3j) 

d($I) = 

= {d t 4> + P<P + vV x <p + uA{D)<f)} lxM dt + y/2uI((di<f>)a^x(a,t)dW p . 



6 



(32) 



Using (|12|) and (|14|) we have 

= I {-2is(d k (a kp ))a jp (d 3 (f)) - (div x U)(p} lx{a)t) dt+ 
+V2uI({dj(f>)ajp)ix{a,t)dW p . 

Now, by Ito, 

dM = Dd($I) + ®ID + d(D, $7). 
In view of (119)1 and ()32J) we have 

d(D,M) = 2i/ J DJ{(a fc cr fep )a- jp (a 7 -0)} |XM) dt 

and consequently the terms i:2uDI(d k a kp )aj P (dj<p)dt cancel and we obtain 

dM = ID {-(div w U)(j) + (div x v + 2vE)}, x , t) dt+ 
+V2uID {a jp (dj(f)) + {d k ^kp)} lX (a,t) dW v 

Now, in view of ((H)) we have that 



j pi 



(div x v) - (div x U) = vdj[{d k a jp )(j kp } - udj[(d k a kp )a. 
and therefore the coefficient of dt in dM is 

DI{2vE + vdj[{d k a jp )(T kp ] - vdj[(d k a kp )(Tj P ]} 

Now 

DI {vdj[{d k a jp )a kp } - vdj[(d k cr kp )a jp ]} = 
DI {p(d k a jp )(d j a kp ) - v{d k a kp ){d j a jp )} = 
= Dn J2 k <j E P Mdk(r jp )(dja kp ) - v(d k a kp )(djCT jp )} = 

= -2vE 

and therefore the coefficient of dt in dM vanishes. We obtained 

dM = V^ID {a ]P {8 3 (j)) + {d k a kp )} {x(a t) dW p , (33) 
that is, M is the martingale 

M(a,t) = 

0(a, 0) + V^i/JjJ /(a, s)D(a, s) {(r jp (dj(f)) + (d k cr kp )} lx(a s) dW p (s) 



7 



Theorem 4 Let h and po be smooth time independent deterministic func- 
tions. Consider the stochastically passive scalar 9 ho (x,t) = h (A(x,t)) and 
the process ip Po of \2b\) with initial datum po- Consider also <f)(x,t), a deter- 
ministic solution of Alty) . Then the random variable 



£{t) = / (f>(x,t)ip po (x,t)9 ho (x,t)dx (34) 
is a martingale. In particular 

E(£(t))= / (j)(a,0)p (a)h (a)da (35) 

holds. 

Proof. In view of the change of variables formula and the definition of ifj po 
we have that 

£ (t) = / M(a,t)po(a)h (a)da (36) 

with M given in (J31|) . The result follows then from the previous proposition. 
More precisely 

d£ = 

V2u{J Rn exp {f° P(X(a, T),r)dT\ a=A {x,s)} Wjpidjfy + (d k a kp )} dx) dW p 

(37) 

gives explicitly the SDE obeyed by £. 

5 Generalized Relative Entropies 

We take now a smooth deterministic, time independent function H of one 
variable, a deterministic solution of (fT2j) . two smooth deterministic, time 
independent functions fo and po, of which p is strictly positive. We form 
the processes ijj po and ipf given by the expressions (J2HJ). Then it t follows 
that 
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holds. Thus, the quantity of interest, ijj po (f)H is the product of a 

stochastically passive scalar, ifj po and 4>. By the previous theorem we have 
that 

8{t) = J i> P0 (x,t)H (|^|) 4>{x,t)dx (38) 
is a martingale. The expected value is then constant in time: 

d Il\ [i> po H(p-)<j>dx\=0. (39) 



dt [J VVw j 

If we denote 

f(x,t)=Ei(> fo (x,t) (40) 

and 

p(x,t)=Ei{> po (x,t) (41) 

we have from (J27)) that / solves (jTHjl . p > solves (jTTjl . We prove that we 
have lfI3|) . 

The starting point is ([3*9^1. In view of (j4T)j) and (jHJ), the statement that 
needs to be proved is 

/ E(W " (Hi) 0dx £ E (S W (42) 

The conservation (J39*|) works for any if, but we expect (|4"2|) to hold only for 
convex if. Indeed, (jHJ) can be reduced to a Jensen inequality. We claim 
more, that for all x, t we have 

Considering the functions 



W P0 ) 

and 



we see that (I43j) becomes 



(44) 
(45) 



H(E(v))<E{gH[- ) }. (46) 
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This, however, is nothing but Jensen's inequality for the probability measure 

Ph = E(gh), 



H P 



< PH 



A Proof of Theorem El 

We devote this appendix to proving Theorem El The original proof can be 
found in [3 j for constant coefficients, and in jl] for variable coefficients. 

Lemma 1 Let X be the stochastic flow defined by (jTSJl . (fTBjl. Then the map 
X is spatially smooth (almost surely), and the determinant D = det(VX) 
satisfies the equation 



dD = D 



( V • v + 2uE) dt + \/2vd k a kv dW } 



where 



E = \ [diO-ipdjGjp — djCTipdiO-jp} . 
Proof. Differentiating (j!5|) we have 

d(d a Xj) = d k Vjd a X k dt + V2ud k a jp d a X k dW p . 



(47) 



Let S n be the permutation group on n symbols, and e T denote the signature 
of the permutation r G S n . By Ito's formula, 

dD= 2 e - Hd c X Tc d(d b X n )+Yl II d d X Td d(d b X n ,d c X Tc ) 

T eS n L c^b c<b d^b,c 

b=l...n 

= e r d k v Tb d b X k Y[ d c X Tc dt + V2vd k a Tb , p d b X k JJ d c X Tc dW p + 

c^b c+b 



T£S n 
b=l...n 



+ v $Z d bXidia TbtP d c X m d m a TcjP J\ 9 d X Td dt 

Cy^b dj^b,c 



(48) 
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We compute each of the terms above individually: 

^2 e T d k v Tb d b X k Y[d c X Tc = ^ trd Tb v Tb d b X Tb JJ d c X Tc + 

reS n c^b reS n c^b 

6=1. ..n 6=1. ..n 

+ ^2 ^e T d k v Tb d b X k d T -iX k [ [ d c X Tc 

reS n k^r b c^r- 1 
6=1. ,.n R 

= (V • v) det(VX)+ 
+ e rd Tk v Tb d b X Tk d k X Tk Yl d c X Tc 

6=l...nre5 n c^b,k 
k+b 

= (V • v) det(VX) + 

The second term above is zero because replacing replacing r with r o {b k) 
in the inner sum produces a negative sign. 
Similarly we have 

^2 e T d k a Tb ^ p d b X k Y[ d c X Tc dW p = d k o Kp det(VX) dW p . 

reS n c+b 
6=1... n 



For the last term in (jlSjl , the only difference is that we have a few extra cases 
to consider: When I = r(6), m = r(c), we will get det(V X)di<Ti p dj(Tjp. When 
I = r(c) and m = r(b), we will get — det(V X)dja ip diaj p . In all other cases 
we get 0. This concludes proof of Lemma [TJ 

Lemma 2 For any time t, the map X t has a (spatially) smooth inverse. 
Proof. Define A by 



A 



exp f (V-v + 2vE - u(d k o- kp ) 2 ) alt + v 7 ^ / d k a kp dWjfi 
Jo Jo 



The Ito's formula immediately shows that A satisfies equation p9jl . Since 
is a linear SDE with smooth coefficients, uniqueness of the solution 
guarantees D = exp (A) almost surely, and hence D > almost surely. 

The spatial invertibility of X now follows as X t is locally orientation 
preserving and has degree 1 (because X t is properly homotopic to X , the 
identity map). The (spatial) smoothness of the inverse is guaranteed by the 
inverse function theorem. 
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The above lemma shows existence of a spatial inverse of X. As before, we 
let A denote the spatial inverse of X. We now derive a stochastic evolution 
equation of A [equation (|2T|)]. 

Lemma 3 Let Y be a C 1 stochastic flow of semi-martingales adapted to Tt, 
the filtration of Wt- If for all a G W 1 , t > we have 

Y(X s (a),ds)= / b(X s (a),s)ds+ / ^.(o),^, 



o 



then 

Y t (a) =Y (a)+ I b(a,t)dt+ I a'(a,t)dW t . 



Proof. Let Y' be the process defined by 

Yj(a)=Y (a)+ [ b(a,t)dt+ [ a\a,t)dW t , 
Jo Jo 

and set 5 = Y — Y' . Since 5 is adapted to !F t , there exists a non-negative 
predictable function a such that 

t 

a(x,y,s) ds= (8(x), 8(y)) t . 



o 



Now, by definition of the generalized Ito integral we have Jq5(X s ,s) = 
almost surely, and hence f Q a(X s , X s , s) ds = almost surely. Since X is 
a flow of homeomorphisms (diffeomorphisms actually), we must have V£, 
a(x, x, t) = almost surely. Thus 5 = Y — Y' is of bounded variation. 
Since we have shown above that 5 has bounded variation, 



/ 8(X s ,ds) = / d t 5\ x s ds 
Jo Jo 



and hence V£, d t S t = 0. At time 0, 5o = by definition, and hence 8 t = 
almost surely for all t, concluding the proof. 

Lemma 4 There exists a process B of bounded variation such that 

A t = B t - V2u [ {VA s )a dW s (49) 
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Proof. Applying the generalized ltd formula to A o X we have 

0= [ t A(X s ,ds)+ fvA\ x s dX s + \ t d%A\ x d{X^,X^) s + 
Jt' Jt' s ' Jt' s ' 

+ / f d i A{X s ,ds),X?-Xf 



A(X s ,ds) + 



VA\ x v + va^ 3 A\ 



ds+ 



(50) 



+ V2u I VA\ a dW s + ( I diA{X s , ds),X? ] - xf 



Notice that the second and fourth terms on the right are of bounded variation. 
Applying Lemma El we conclude the proof. 

Lemma 5 The process A satisfies the equation 



dA t +(vV)A t dt-ua ij df j A t dt-V2ijd j A t (d i a jk )a ik dt+(VA t )adW t = (51) 

Proof. Since the joint quadratic variation term in (|5Up depends only on the 
martingale part of diA, we can compute it explicitly by 

diA(X s ,ds), = -2v (d?Acrjfc0ifc + d j A s (d i a jk )a ik ) oX s ds 

= -2v / (ciijdfjAs + d j A s {d i a jk )a ik ) o X s ds. 
Jt' 

(52) 



Substituting (JH2j) in ((SD| and applying Lemma El we conclude the proof. 
Acknowledgment. PC partially supported by NSF grant DMS-0504213. 
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